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ABSTRACT 

In a previous paper (Brunetti et al. 2004) we presented the first self-consistent calcu- 
lations of the time-dependent coupled equations for the electrons, hadrons and Alfvcn 
waves in the intracluster medium, which describe the stochastic acceleration of the 
charged particles and the corresponding spectral modification of the waves. Under 
viable assumptions, this system of mutually interacting components was shown to ac- 
curately describe several observational findings related to the radio halos in clusters 
of galaxies. 

In this paper, we add to the self-consistency of the calculations by including the 
generation and re-energization of secondary electrons and positrons, produced by the 
inelastic interactions of cosmic rays with the thermal gas in the intracluster medium. 
The bulk of Cosmic rays is expected to be confined within the cluster volume for 
cosmological times, so that the rate of production of secondary electrons, as well as 
gamma rays, may become correspondingly enhanced. If MHD waves are present, as 
it may be expected in the case of a recent merger event, then the reacceleration of 
secondary electrons and positrons can significantly affect the phenomenology of the 
non thermal processes in clusters. We investigate here these effects for the first time. 

Key words: acceleration of particles - turbulence - radiation mechanisms: non- 
thermal - galaxies: clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

The intracluster medium contains a non thermal compo- 
nent in the form of magnetic fields and relativistic particles, 
as demonstrated mainly by the appearance of diffuse radio 
emission (e.g., Feretti, 2003) and by studies of the rotation 
measure of radio sources in galaxy clusters (e.g., Clarke et 
al., 2001), but also by the hard X-ray (HXR) excess de- 
tected in a few galaxy clusters by the BeppoSAX and RXTE 
satellites (Fusco-Femiano et al. 2003; Rephaeli & Gruber 
2003). While the radio radiation is certainly the result of 
synchrotron emission of high energy electrons in the intra- 
cluster magnetic field, HXRs may be explained in terms of 
inverse Compton scattering (ICS) of the same electrons ofT 
the photons of the cosmic microwave background (CMB), or 
as a result of bremsstrahlung emission from supra-thermal 
electrons (e.g., Ensslin, Lieu, Biermann 1999; Blasi 2000; 
Dogiel 2000; Sarazin & Kempner 2000). The latter solution 
does however require an energy input which may be sup- 
ported only for a time shorter than a few hundred million 



years, otherwise the gas is overheated (Petrosian 2001; Blasi 
2000). It should be mentioned that the rather poor sensitiv- 
ity of the present and past facilities for the observation of 
HXRs does not allow to obtain a iron-clad detection of HXR 
excesses. Future observatories (e.g. ASTR0-E2, NEXT) are 
necessary to confirm (or disprove) the existence of these 
excesses (see Rossetti & Molendi 2004; Fusco-Femiano et 
al. 2004). 

So far the most serious problem from the theoretical 
point of view is to identify the origin of the radiating high 
energy electrons. Buote (2001) and Schuecker et al. (2001) 
have found hints of a correlation between the non thermal 
diffuse radio emission and the presence of merger activity 
in the host clusters. This may indicate a link between the 
process of formation of galaxy clusters and the origin of the 
non thermal activity. 

Two main avenues have been identified to explain the 
fact that such high energy electrons are present and able 
to radiate on distance scales larger than their typical loss 
lengths: in the context of the so-called primary models elec- 
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trons are accelerated at shock waves through the first order 
Fermi mechanism or they are continuously re-energized vn 
situ on their way out (Jaffe 1977); the latter case is also 
called reacceleration model. On the other hand, in the con- 
text of secondary electron models electrons arc secondary 
products of the hadronic interactions of cosmic rays with the 
intracluster medium, as first proposed by Dennison (1980). 

Shock waves are unavoidably formed during merger 
events due to the supersonic relative motion of two (or more) 
infalling clusters. Since electrons have a short pathlength 
due to IC losses, they can only move a short distance away 
from the acceleration regions. If this acceleration site were 
to be at a shock front, the emission from the electrons would 
be concentrated around the shock rim (e.g., Miniati et al. 
2001) and the spectrum of the radiation would be quite steep 
because of the low Mach numbers in the central virialized 
regions of clusters (Gabici and Blasi 2003, Berrington and 
Dormer 2003), where the radio emission is observed to come 
from. On the other hand, strong shocks arc formed in the 
outskirts of clusters, and they may be responsible for the 
acceleration of protons (Ryu et al. 2003). These protons can 
then be advected into the cluster and there be confined for 
cosmological time scales (Volk et al. 1996; Berczinsky, Blasi 
& Ptuskin 1997). The confinement unavoidably increases the 
energy density of cosmic rays in the intracluster medium, 
and may correspondingly favour the generation of gamma 
rays through the decay of neutral pions and of electrons and 
positrons through the decay of charged pions. 

The origin of the emitting particles is still matter of 
debate (e.g., Ensslin 2004). The different models for the 
production of the radiating electrons have however a sub- 
stantial predictive power, which can be used to discriminate 
among such models by comparing their predictions with ob- 
servations. Recently Brunetti (2004) and Blasi (2004) have 
discussed at length the strong and the weak points of the 
models in explaining the phenomenology of existing data. 
Although the role of future observations remains crucial in 
order to achieve a definite conclusion of the right descrip- 
tions of the non thermal phenomena we detect, both au- 
thors concluded that, at least as far as the Coma cluster 
and few other well studied clusters are concerned, present 
data seem to suggest the presence of particle-reacceleration 
mechanisms active in the ICM. On the other hand, several 
complex and poorly understood processes are involved in 
these models. 

Cluster mergers induce large-scale bulk fiows with ve- 
locities ~ 1000 km s~^ or larger. These flows drive insta- 
bilities on large-scales which redistribute the energy of the 
mergers through the cluster volume and decay into turbu- 
lent velocity fields. It has been shown that reacceleration 
of a population of relic electrons by turbulence powered by 
cluster mergers is a promising mechanism to explain the 
very large scale of the observed radio emission, the complex 
spectral behaviour observed in some diffuse radio sources 
(Brunetti et al., 2001; Petrosian 2001; Ohno, Takizawa and 
Shibata 2002; Pujita, Takizawa and Sarazin 2003), and the 
observed occurrence of radio halos with cluster mass (Cas- 
sano & Brunetti 2005). 

A step forward in the study of particle acceleration in 
gala^xy clusters has been achieved by recent studies of the 
interaction between particles and Alfven waves in a very 
general situation in which relativistic electrons, thermal pro- 



tons and relativistic protons are present in the cluster vol- 
ume (Brunetti et al., 2004, hereafter Paper I). In Paper I 
the interaction of all these components with the waves, as 
well as the turbulent cascading and damping processes of 
Alfven waves, have been treated in a fully time-dependent 
way in order to calculate the spectra of electrons, protons 
and waves at any fixed time. This represented the first at- 
tempt to include in a self-consistent way hadronic cosmic 
rays in the reacceleration scenario. It was found there that 
radio halos and HXR tails could be activated for a time scale 
of ~ 0.5 — 1 Gyr through resonant interaction of particles 
with short-wavelength Alfven waves, coming from the de- 
cay of merger induced turbulence. The role of the liadrons 
was found to be that of exerting a substantial backreaction, 
so that the non thermal activity is suppressed if the en- 
ergy content in the form of hadrons gets larger than a few 
percent of the thermal energy. We named this phenomenon 
wave-proton boiler. 

The present paper serves as a completion of the effort 
started in Paper I: here we add to the previous calcula- 
tion the effect of the reacceleration of those electrons and 
positrons that are generated at any time as secondary prod- 
ucts of the inelastic interactions of the confined cosmic rays 
with the thermal gas. 

In Sections 2 and 3 we provide the reader with a review 
of the main aspects of the calculations developed in Paper 
I. In Section 4 we describe the general scenario in which sec- 
ondary reaccelerated electrons arc present. The application 
of the full calculation to the phenomenology of a Coma-like 
cluster is presented in Section 5. We conclude in Section 6. 



2 PARTICLE ACCELERATION AND ENERGY 
LOSSES IN THE ICM 

In this Section we briefly summarize the rates of energy loss 
and gain for non thermal leptons and protons in the intrar 
cluster medium, relying on the formalism already introduced 
in Paper I. 



2.1 Energy losses for leptons 

Four channels dominate the energy losses of relativistic lep- 
tons with momentum p, namely ionization. Coulomb scat- 
tering, synchrotron emission and IC. The rate of losses due 
to the combination of ionization and Coulomb scattering can 
be written following (Sarazin 1999): 



\dtJi 



-3.3 X 10" 



nth 



1 + 



ln(7/nth) 
75 



(1) 



where nth is the number density of the thermal plasma. The 
rate of synchrotron plus IC losses is: 



\dtj rad 



4.8 X 10"V 



(If)' 



2/3 



(2) 



where B^g is the magnetic field strength in units of jj,G, 
and G is the pitch angle of the emitting leptons; in case of 
efficient isotropization of the electron momenta, the sin^ 9 
is averaged to 2/3. In the typical conditions of the ICM, 
radiative losses are dominant for leptons with Lorentz factor 
7 >> 100, while Coulomb losses dominate at lower energies 
(Sarazin 1999,2002; Brunetti 2003). 
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2.2 Energy losses for protons 

For relativistic protons, the main channel of energy losses 
in the ICM is provided by inelastic proton-proton collisions. 
The time-scale associated with this process is: 



1 



nthfppC 



10' 



i( nth \ ^ 



(3) 



Thus inelastic pp scattering is weak enough to allow the ac- 
cumulation of protons over cosmological times (Borozinsky, 
Blasi and Ptuskin 1997), however it is also efhcicnt enough 
for the continuous production of pions, which in turn decay 
into gamma rays (for neutral pions) , electron-positron pairs 
and neutrinos (for charged pions). The process of pion pro- 
duction in pp scattering is a threshold reaction that requires 
protons with kinetic energies larger than about 300 MeV. 

For trans-relativistic and mildly relativistic protons, 
energy losses are dominated by ionization and Coulomb 
scattering. Protons more energetics than the thermal elec- 
trons, namely with 0p > Be = (3/2me/mp)"'"''^/3e (,5e — 
0.18(T/10*/s:)^/^ is the velocity of the thermal electrons), 
arc affected by Coulomb interactions. Defining Xm = 



(^) 



1/3 



/3e, one has (Schlickeiser, 2002): 



f ^-1.7x10-^^ 
at 



f nth \ 
UO-sJ 



/3p 



cgs units, 



which has the following asymptotic behaviour 




mcPc < p < mcxm 
mcxm < p « mc 
p » mc 



(4) 



(5) 



2.3 Alfvenic acceleration of relativistic particles 

Different coupling between particles and waves may result 
in energy transfer from magnetic fluctuations into relativis- 
tic particles: Magneto-Sonic (MS) waves, magnetic Landau 
damping (Kulsrud & Ferrari 1971; Schlickeiser & Miller 

1998) , Lower Hybrid (LH) waves (e.g., Eilek & Weatherall 

1999) and Alfvcn waves arc a few examples that have been 
investigated in the literature. Alfven waves efficiently couple 
with relativistic particles via resonant interaction and they 
are likely to transfer most of their energy directly into these 
relativistic particles. 

The resonant condition for a wave of frequency ui and 
wavenumbcr projected along the magnetic field fcy, for a 
particle of species a with energy Ea and projected velocity 
uy = vn is (Melrose 1968; Eilek 1979): 



(6) 



1 



where, in the quasi parallel case {k± « ma^ia/p), v 
(y = 1) for electrons (protons and positrons). 

The dispersion relation for Alfven waves in an isotropic 
plasma with both thermal and relativistic particles was 
given by Barnes & Scargle (1973). In the conditions typical 
of galaxy clusters, the dispersion relation of Alfven waves 
reduces to ~ |fc|| |ua. Combining the dispersion relation of 
the waves with the resonant condition, Eq. 6, one can de- 
rive the resonant wavenumber, kres, for a given momentum, 
p = mvy, and angle, //, of the particles : 



(7) 



where the upper and lower signs refer to protons (and 
positrons) and electrons respectively. The interaction of par- 
ticles with Alfven waves can be thought of as a diffusion 
process in the momentum space of the particles. If the dis- 
tributions of waves and particles are assumed to be isotropic, 
then the diffusion coefficient was found by Eilek & Henriksen 
(1984) to be: 



Dpp{p,t) 



o 2 2 2 

271 eu^ 



Wkjt) 
k 



\ c ^ pk J 



dk,{8) 



where the minimum wavenumber (maximum scale length) of 
the waves interacting with particles with given momentum 
is: 

Om 1 
^ M± 

and kmax is given by the largest wavenumber of the Alfven 
waves, which is fixed by the condition that the frequency 
of the waves cannot exceed the proton cyclotron frequency, 
namely a; < Sip. It follows that kmax ~ 0.p/vA (in the follow- 
ing, for consistency with Paper I, we take kmax ~ ^Ip/vM, 
vm being the magnetosonic velocity). 



3 EQUATIONS AND COUPLING 

3.1 Time Evolution of Particles and Waves 

In this Section we summarize the formalism that we intro- 
duced in Paper 1 for the description of the time-dependent 
interaction between the particles and the waves. The evolu- 
tion of the electron (and positron) number density is given 
by the diffusion equation, which includes energy losses and 
gains (Paper I and references therein): 



dN,{p,t) ^ d_ 
dt dp 



Ne{p,t) 



dp 

dt rad 



dp 
'dti 



p-' op 



dp^ 



+ [DppN,(p,t)\ + Qe\p,t;Np{p,t)\. (10) 



Here Qe[p,t; Np{p,t)] represents the injection rate of sec- 
ondary relativistic electrons and positrons generated dur- 
ing the collisions between the accelerated relativistic protons 
with the thermal protons in the ICM (Sect. 3.3). A similar 
equation can be written for protons (Paper I and references 
therein) : 



dNp(p,t) ^ d_ 
dt ~ dp 



Np{p,t) 



dti 



=l2 



dp- 



■[DppNp(p,t)]. (11) 



The evolution of the spectrum of the Alfven waves is 
described through a diffusion equation in the wavenumber 
space (Eilek 1979): 



dWkjt) _ _a / dW^jt) 

dt '^'^ dk 



J2^\k)W^{t) 
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(12) 



where -Dkk is the diffusion coefficient due to wave-wave cou- 
piing (Sect. 3.2.2), F is the damping rate of the Alfvcn waves 
with particies (Sect. 3.2.1) and Ik is the injection rate of the 
Alfven waves (Sect. 3.2.3). In Eq.(12) we use the assump- 
tion, commonly made, that wave-wave interaction is just 
local in the wave number space. 



3.2 Damping Procesess, wave-wave coupling and 
injection of turbulence 

All the relevant processes related to wave-particle interac- 
tions, wave-wave interactions and injection of Alfven waves 
in the ICM are described in detail in Paper I. Here we only 
provide the reader with a brief overview of the main pro- 
cesses that we include in our calculations. 



3.2.1 Damping Processes 

In the case of nearly parallel wave propagation (i.e., k± « 
mQ/p, ~ feji) and isotropic distribution of the velocities 
of the particles, the cyclotron resonant damping rates for 
Alfven waves with particles of species a are given by Melrose 
(1968): 



2/1 .2^dfc{p,t) 



dp 



dp ■■ 



2 2 2 

n e Va 



{i-^l)h^^^-'J^\dp, (13) 



kc^ .L ^ \ P dp 

where, in the relativistic case one has : 

rel VA , fla^a 
Ha = ± 

c pk 
and in the non relativistic case: 

th VAaia , ^laBla 
Mo = ± 



P 



pk 



(14) 



(15) 



The upper and lower signs in Eqs. 14-15 are for negative 
and positive charged particles respectively. Eq. 13 can also 
be used to evaluate the damping rate in the case of isotropic 
Alfven waves with an approximation which is within a factor 
of ~ 3 (Lacombe 1977). 



3.2.2 Wave-wave Cascade 

Wave-wave interactions cause the spectrum of the waves to 
cascade, namely to broaden toward larger values of k. This 
is a diffusive process, with diffusion coefficient -Dkk = k^ /tb. 
The time Ts is the spectral energy transfer time and can be 
written as Tg ~ t^i^/ts (Zhou & Matthaeus 1990), where 
tml = \/Sv is the non- linear eddy-turnover time {5v is the 
rms velocity fluctuation at wavelength A) and ts is the time 
over which this fluctuation interacts with other fluctuations 
of similar size. 

In the context of the Kolmogorov phenomenology, the 
Alfven crossing time ta ~ X/va exceeds tnl and the fluc- 
tuations of comparable size interact in one turnover time, 
namely ts ~ tnl. Since the velocity fluctuation, 5v, is re- 
lated to the rms wave field, SB, by Sv^/va = 5B^/B^, the 



diffusion coefficient can be written as (Miller & Roberts 
1995): 



Dkk ~ va/c 



7/2 



WkW 

2Wb 



1/2 



(16) 



Given a spectrum of injection of waves per unit time, Ik, 
one simple possibility to estimate the cascade time scale is 
to use the spectum of the waves in Eq. 16 as obtained from 
Eq. 12 under stationary conditions and without damping 
processes. In Paper I we found: 



Dkk 



1 /B^ 
fc V 47r 



(17) 



3.2.3 Injection of Alfven waves in the ICM: the Lighthill 
mechanism 

While the physics involved in the process of energy trans- 
fer between waves and particles for a given spectrum of 
waves is relatively well understood, the transformation of the 
wave spectrum starting from some injection at large scales 
is rather poorly known. The waves are expected to couple 
with relativistic particles when the turbulence has been en- 
riched of short wavelength modes, so that the cascading is 
implicitely required to be rather eflScient if the injection oc- 
curs on macroscopic scales. If however this is the case, it 
was shown (Yan & Lazarian 2004 and refs. therein) that the 
Alfven waves reach the high-k part of the spectrum with 
a highly anisotropic spectrum, and the efficiency of parti- 
cle acceleration is likely to be therefore drastically reduced. 
From this follows that the acceleration process is favored in 
those scenarios in which the injection of Alfven waves oc- 
curs on relatively small scales to start with. One injection 
process in which this condition is fulfilled is provided by 
the so-called Lighthill mechanism (Kato 1968; Eilek & Hen- 
riksen 1984), which can convert some fraction of the large 
scale fluid turbulence on the larger scales into Alfven waves 
on smaller scales. Following Fujita et al. (2003) and Paper I, 
we assume in our calculation that fluid turbulence is injected 
on large scales, for instance excited by a merger event, and 
that the Lighthill mechanism couples the fluid turbulence 
with MHD turbulence on smaller scales. We made the as- 
sumption here that the spectrum of the fluid turbulence (not 
the MHD turbulence) is in the form of a power law 



W({x{) = WiXj' 



(18) 



in the range Xo < xi < a;™'^, where Xo is the wavenum- 
ber corresponding to the maximum scale of injection of the 
turbulence; we thus do not consider the possibility that tur- 
bulence may be injected simultaneously at many scales (e.g., 
Tsytovich 1972). The maximum wavenumber is where the 
effect of fluid viscosity starts to be important and it is of 
the order of a;^^'' ~ x^iTlf^'^ (Landau & Lifshitz, 1959), 
TZ being the Reynolds' number. Thus for high values of the 
Reynolds number in the ICM (and in a magnetized medium) 
the turbulence cascading can be an efficient process down to 
scales much smaller than kpc (Fujita et al. 2003; Paper I). 

In the Lighthill process a fluid eddy may be thought 
of as radiating Alfven waves at a wavermmber k = 
(vi(x)/vA)xi. The Alfven waves are expected to be driven 
only for a;/ > XT, xt being the wavenumber at which the 
transition from large-scale ordered turbulence to small-scale 
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disordered turbulence occurs. This transition is usually as- 
sumed to take place at the Taylor scale (Eilek & Henrik- 
sen 1984), It ~ lo{15/7iy^^ , where the Reynolds number is 
given by 7?. = loVf /i/k , and i>k is the kinetic viscosity. 

The fraction of the fluid turbulence radiated in the form 
of MHD modes is small for all but the larger eddies, near 
the Taylor scale. Therefore the Lighthill radiation may be 
expected to not disrupt the fluid spectrum. The rate of rar 
diation via the Lighthill mechanism into Alfven waves of 
wavcrmrnbor k is (Eilek & Henriksen 1984; Pujita et al.2003; 
Paper I): 



Ik 



2m 



3 

PVA 



\vlRj 



X k 



(19) 



where p ~ St/v'f, with £t the energy density of the fluid 
turbulence, and 



R = 



XoWf{Xo) 
XTWf{XT) 



(20) 



3.3 Secondary Electrons 

As discussed above, the main new ingredient added in this 
paper, compared with the calculations presented in Paper 
I is the presence of secondary electrons (and positrons), as 
generated in the hadronic inelastic interactions of cosmic 
rays with the thermal gas in the ICM. The decay chain that 
we consider is (Blasi & Colafrancesco 1999): 



P + P 



n" + Tv^ + TV + anything 



77 



The spectrum of secondary electrons and positrons with en- 
ergy Ee is given by the convolution of the spectra of pro- 
tons, N{Ep), with the spectrum of pions produced in a single 
cosmic ray interaction at energy Ep, (FTr{ETr,Ep)) and with 
the distribution of leptons from the pion decay, F^{Ee, -Btt), 
(e.g., Moskalenko & Strong, 1998): 

Qf\p,t;Np] = vF^f^c [ dEpf3^N{Ep)af{Ep) 



dE^F^{E^, Ep)Ft{Ee, E„), 



(21) 



where a^{Ep) is the inclusive cross section for pion produc- 
tion, Etr is the threshold energy for the process to occur and 
the distribution of electrons and positrons is given by : 



Ft{Ee,E^±) — j dEfjtF^ {Ee, Efj,, Eir), 



(22) 



where F^{Ee, E^, Ett) is the spec- 

trum of electrons/positrons from the decay of a muon of 
energy produced in the decay of a pion with energy S„ . 

At large values of Ep the differential cross section is 
sufficiently well described by the so-called Fcynman scaling, 
with small deviations which can easily be taken into account. 
In the low energy part, when the reaction occurs close to 
the threshold, and in general at laboratory energies smaller 



than ~ 10 GeV, the experimental data on pion produc- 
tion arc rather poor, and the scaling behaviour is violated. 
Since in this paper we axe going to calculate the spectrum 
of the reaccelerated secondary electrons and positrons, we 
arc forced to use a source term which correctly describes the 
spectrum of the injected leptons over a broad energy range 
(7 ~ 10^ — 10^). A practical and useful approach to both the 
high energy and low energy regimes was proposed in Dermer 
(1986a) and reviewed by Moskalenko & Strong (1998), and 
is based on the combination of the isobaric model (Stecker 
1970) and scaling model (Badhwax et al., 1977; Stephens 
& Badhwar 1981). Here we briefly describe the formalism 
and approximations used in our calculations and provide 
the main equations. 

In the Stecker's model the pion production due to pp 
collisions near threshold is mediated by the excitation of 
the A3/2 isobar, which subsequantly decays into a nucleon 
and a pion. In this case the spectrum of pions produced in a 
single cosmic ray interaction is given by (e.g., Dermer 1986a; 
Strong & Moskalenko 1998): 



^7r(-E'x, Ep) — P 



^ """" dmAC^U{E^,Ep-mA) 
2^ 2 (mAc2-m°c2)2-KP2 

— rUpC' 



"2 -m^c^ 



tan-i { m,c'+m^f-m%^ \ , ^^g^ 



where P ~ 0.0575 GeV is the width of the Breit-Wigner 
distribution, rn^c^ ~ 1.232 GeV is the average rest energy 
of the A-isobar, and s = 2mpC^{Ep + mp(?) is the square of 
the energy in the center of mass frame, 



fn{E^,Ep;mA) 



H+ H- 
+ 



(24) 



where = 1 for 7±7;(1 - /3±/?;) < 7, < 7^(1 + /3±/3;) 
and = otherwise, 

7a = 7c7l(l ± PcPl), 



with 



7a = 



, 2 4 2 4 

s + m^c — rUpC 

2y^mAc2 



(25) 



(26) 



the Lorentz factor of the A-isobar in the center of mass 
frame, and 



7c = y/sj^mpC 

is the Lorentz factor of the center of mass. Finally, 



7vr = 



mA + fi^ — rr?p 



(27) 



(28) 



2 m A w^TT 

is the Lorentz factor of the pion in the A-isobar system. 

At high energies the spectrum of pions can be approx- 
imated by the simple formula proposed by Berezinsky & 
Kudryavtsev (1990): 



Fn{ETr, Ep) — - 



ci(l-|l) +C2exp(-18|l) 



, (29) 



where ci = 1.22 and C2 = 0.92. Thus the injection rate of 
pions is given by : 
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and the injection rate of relativistic electrons/positrons is 
given by : 



F.^{E^,Ep)a^{pp 



■,(30) 



for f(i+/3.r< 2^^^^ 

Here = 2m'iEe/{m'i + m^)S,r, and we put 



^^ = -7^^^{4-[l+(^) Tj- (37) 
r Pit ml - mt L V rrt^ / J J 

Qe±{P.t)= I QAK±,t)dE^ dEi,F^±{E^,Ei„E^)Fi,{E^,E^),{31) 

J Eqs. (36-37) are then combined with Eq. 10 to calculate the 



where, following Moskalenko & Strong (1998), in the calcula- 
tion of the pion injection rate we combine the isobaric model 
(Eqs. 23-24) with the scaling model (Eq. 29) and adopt a 
linear interpolation between the two regimes, in the energy 
range 3-7 GcV. In our calculations wo adopt the fits to the in- 
clusive cross section a^{Ep) given in Dermer (1986b) which 
allow to describe separately the rates of generation of tt"*" 
and 7r~. 

The pion decay is well known to generate a muon spec- 
trum in the following form: 

™2 



Fi^{Ei^,ETr) - 



1 



ml -mj^pn' 



(32) 



time evolution of the spectrum of the accelerated secondary 
leptons. 

Finally, it is useful to derive the injection rate of sec- 
ondary electrons/positrons at high energies (in the scaling 
approximation) and for a simple power law spectrum of cos- 
mic ray protons, N{Ep) = KpE~^. In this case the inclusive 
cross section for the production of tt^ is approximately that 
of 7r° and thus from Eqs. (36-37), one finds : 

QMp)=Ais)n,,KpE:^[p.+V.ln(^)+Vs (|^) ^ } 

where A{s) = 64 x c x 10~^'^a^~^, and we put 

2 ml 



Muons are produced in a relatively narrow range of energies, 
between a kinematic mimimum and maximum given by 



(39) 



-, Ett 

and 



2mi 



(m^(l-/3^)+m^(l + /3.)), 



(m^(l+/?.)+m^(l-/3.)) 



(33) 



(34) 



and 
Pi = 

V2 = 



+ 



(s + 2)2 (s + 3)^ 



To 



h 



In order to speed up the computation, we assume that the 
spectrum of muons is a delta-function at the energy E^ = 
l/2(£j^,niin -I- -B^,max). Therefore: 



(35) 



Vz = 1.5 X 



s4-2 sH-3 



-I- 



F^iE^,E^)^5\ E^-E. 
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where 
5 

a = 



The spectrum of electrons and positrons from the muon 
decay, F^±{ETr,E^,Ee), was calculated by Blasi & Co- 
lafrancesco (1999). Combining their results with Eqs. (31) 
and (35), we obtain the rate of production of secondary elec- 
trons/positrons : 



s + 1/2 s + 5/2 s + 7/2 



12 

6=-|(l + (a^-l) 



Qe±iP,t) 



where 



1 + 2{a - 1] 



1 + (mp/m^)^ 
1 - (m^/m^)2 

1 + (m^/m^)^ 
1 — {m^/m^y^ 



dEiT ^ „ Fe{Ee, ETr),{S6) and whore the integrals* are defined as: 
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ci(l - x)i + C2 exp(-18a:) 
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s+2 s+3 
(41) 

l2 (42) 

(43) 
(44) 
(45) 

(46) 
(47) 

(48) 



3(1 - /3.) 



MX 



It follows that the slope of the spectum of secondary elec- 
trons and positrons is essentially that of the cosmic ray pro- 
tons in this approximation (e.g., Dermer 1986a; Blasi & Co- 
lafrancesco 1999). In addition, Eq. 38 describes the slight 



* Pi=0.069, 0.046, 0.031, 0.022, ^2=0.019, 0.013, 0.01, 0.007, 
and 7^3=0.006, 0.004, 0.0033, 0.0025, in the case s=2.0, 2.1, 2.2, 
and 2.3, respectively 
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departure from the simple power law shape which is due to 
the increase of the inclusive cross section with the energy of 
the scattering protons (e.g., Dermer 1986b). 

As shown in Paper I, the spectra of protons as affected 
by the reacceleration are usually not power laws. Therefore 
the expressions given here for the case of power law spectra 
are in general not applicable, although in the following we 
will sometimes use them, where specified, in order to esti- 
mate orders of magnitude. 



4 NUMERICAL CALCULATION OF THE 
SPECTRA OF PARTICLES AND WAVES 

4.1 Basic Assumptions 

A detailed modelling of the injection of MHD turbulence 
in galaxy clusters and of all the related processes of wave- 
particle coupling is a very complex matter and well above 
the capabilities of present numerical simulations and semi- 
analytical treatments. On the other hand, the basics of this 
process can be hopefully understood by making use of viable 
assumptions and simplifications: this is the aim of Sect. 4-5. 

We assume that the injection of turbulence starts in 
coincidence with a merger event and remains constant for 
the duration of such an event. As a necessary simplifica- 
tion the spectrum of fluid turbulence is taken in the form 
of a power law (Eq.18), which basically means that there 
is roughly a single driving scale. The turbulence that is in- 
jected is only fluid turbulence, while the MHD turbulence is 
developed later as a consequence of the Lighthill mechanism. 
We assume that the physical conditions in the ICM (namely 
magnetic field strength, temperature and number density of 
the thermal particles) do not change significantly during the 
time in which the turbulence is injected. 

We also assume that there is no spatial diffusion of the 
particles during the period of injection of the fluid turbu- 
lence, and ignore the effect of the mixing processes which 
may take place during cluster merger events. 

In addition we assume that the fluid turbulence and the 
MHD turbulence are isotropic and that the magnetic field 
is tangled enough to ensure that also the distribution of the 
accelerated particles is isotropic in pitch angle. 

With these assumptions the interaction between waves 
and particles can be investigated by solving the set of cou- 
pled differential equations Eqs. 10, 11, and 12. We consider 
situations in which the amount of energy injected in the 
form of turbulence is typically much smaller than the ther- 
mal energy of the ICM, and thus the thermal distributions 
of electrons and protons are treated as stationary. Since the 
time scale of damping and cascading are much shorter than 
the particle acceleration time scale, following Paper I we 
adopt a quasi stationary approach, in which it is assumed 
that within each time-step the spectrum of the waves ap- 
proaches a stationary solution (obtained by solving Eq. 12 
with dW/dt = 0) and that this solution changes with time 
due to the evolution of the spectrum of the accelerated elec- 
trons and protons. 



4.2 The wave-proton Boiler 

As discussed in Paper I, Alfven waves channel most of the 
energy into relativistic protons, therefore subtracting it from 
the electron component which is on the other hand respon- 
sible for the observed radiations. It follows that if protons 
are too abundant in the ICM, the MHD turbulence is too 
efficiently damped and the acceleration of electrons is sup- 
pressed. This process is in fact quite complex, since it de- 
pends on the non linear time-dependent interplay among the 
spectra of protons, electrons and waves. 

The aim of this Section is uniquely to draw a few gen- 
eral conclusions on the efficiency of electron and positron 
acceleration in the ICM. To achieve this goal, we simply as- 
sume that the injection rate of Alfven waves can be written 
as a power law, I{k) = lok^^ , and that, for simplicity, the 
spectrum of relativistic protons can also be approximated as 
a power law, Np{p) oc p"". Neither one of these assumptions 
is adopted in the detailed calculations that follow. 

Within the context imposed by these simple assump- 
tions, the efficiency of lepton acceleration can be estimated 
analytically. In particular, if Ts is the time for spectral energy 
transfer due to wave-cascading and Td is the damping time 
of waves on the relativistic particles, one has two relevant 
asymptotic cases: 

i) For Ts « Td, the spectrum of the waves is driven by the 
injection of Alfven waves and the process of wave-cascading. 
For roughly stationary regime, the spectrum of the waves 
can be estimated from Eqs. (12 and 16), and we can write: 

W, = (nthmp)i (49) 

The lepton diffusion coefficient in momentum space is ob- 
tained by combining Eq. 8 and Eq. 49: 

B'/\ (50) 

nth 

The efficiency of electron and positron acceleration in- 
creases when lo/nth increases and it slightly increases with 
increasing B. Thus, if the turbulence is smoothly injected 
in the cluster volume with an injection rate which scales 
with the thermal energy density, the efficiency of lepton ac- 
celeration is expected to be slightly higher in the high field 
regions. On the other hand, if the injection rate of the tur- 
bulence does not depend on the local energy density of the 
thermal ICM, then electrons are accelerated with higher ef- 
ficiency in the low density regions. 

ii) For Ts >> Td, the spectrum of the waves is determined 
by the injection of Alfven waves and their damping, which is 
mainly caused by resonant scattering with relativistic pro- 
tons. From Eq. 12 and from the damping rate due to relar- 
tivistic protons (Eq. 58 in Paper I) one has: 

W, = ' ( P'°w/(mpc) y- 

The diffusion coefficient in momentum space for electrons 
and positrons can therefore be obtained by combining Eq. 8 
and Eq. 51: 

Dpp oc ^B'-". (52) 

Op 

It follows that the efficiency for lepton acceleration increases 
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Figure 1. Time evolution of the proton spectrum. From bottom 
to top, we plot the spectra for t = (dotted line), 10^®, 3 x 10^^, 
7x 10^^ , and W^^ seconds after the start of the acceleration phase. 
The calculations are carried out assuming n^;, = 1.8x 10~^ cm ~^ , 
T = 10*K, B{r = 0) = 1.4^tG, d/dt{SB)'^/SiT = 1.2 x lO-^^erg 
cm^'^ s^-', and £p = lO^'^ft^ (with s = 2.2). In the Figure we 
also show the axis with the values of the Lorentz factors of the 
secondary electrons and positrons produced by the corresponding 
cosmic ray protons. 



of Paper I, in which wc showed that the spectrum of rela- 
tivistic protons may be substantially modified iu the energy 
range 1 GeV - 100 GeV due to the resonant interaction with 
MHD Alfvenic turbulence. This is a consequence of the fact 
that the acceleration time for protons is minimum at these 
energies (see Fig. 16 in Paper I). Indeed, at smaller scales 
(which resonate with smaller values of the proton energy) 
Ts » Td and the acceleration time should decrease with 
particle's momentum as (from Eqs. 8-9 & 51): 

racc(p)=pV^^oc-|-ocp-'"-H3^ (53) 

while at larger scales (which resonate with more energetic 
protons) Ts « Td and the acceleration time should increase 
with the momentum of the particles (from Eqs. 8-9 & 49) 
as: 

racc(p) ocpi/^ (54) 

In Fig. 1 we also mark the typical regions of the spec- 
trum of the hadrons which approximatively contribute to 
the injection of the secondary electrons and positrons with 
a given Lorentz factor. The consequence of the decrease of 
the efficiency of Alfvcri acceleration with increasing pro- 
ton energy is that only the amount of secondary elec- 
trons/positrons injected at 7 ~ 10^ — 10^ is expected to be 
significantly increased. On the other hand, the injection rate 
of secondary electrons and positrons with 7 ~ 10*, which 
emit the synchrotron radiation at 0.3-1.4 GHz, is not ex- 
pected to be substantially modified (at least not more than 
a factor of 2-3) by the Alfven acceleration process. 



with increasing Io/£p and with decreasing B (at least for 
UJ > 1). 

In general, the damping of waves due to resonant scat- 
tering with protons and the process of cascading of waves 
contribute rougly equally to the modification of the spec- 
trum of the Alfven waves in the ICM, so that a realistic 
situation can be thought of as intermediate between the two 
regimes discussed above, at least at the beginning of the 
reacccloration phase. 

On the other hand, when rcacccleration has started, en- 
ergy is transferred from waves to protons, and the damping 
due to this process is expected to increase with time. This 
effect is more pronounced where the Alfven velocity is larger. 
As shown in Paper I, a general feature of the reacceleration 
picture in the presence of hadrons is that the acceleration 
efHciency is slightly higher in the low density regions so that 
the radio emission due to the synchrotron radiation of elec- 
trons may have a very broad profile as a function of the 
radial distance from the cluster center. 



4.3 The spectrum of relativistic protons 

Due to the confinement phenomenon and the negligible en- 
ergy losses of hadrons in the ICM, the spectrum of the 
protons after the reacceleration phase is left basically un- 
changed. 

In Fig. 1 we illustrate an example of the time evolution 
of the spectrum of protons assuming physical conditions typ- 
ical of the core of a massive cluster. We confirm the finding 



4.4 Reacceleration of secondary electrons and 
positrons 

The main point of this paper is to include the ofi^ect of the 
reacceleration of the secondary electrons and positrons, as 
generated in hadronic interactions of a time-dependent spec- 
trum of protons. This phenomenon has a twofold effect on 
the observable non thermal radiation from a cluster: first, 
the secondary electrons and positrons add to the pool of 
(primary) electrons that can suffer the re-energization due 
to coupling with waves. Second, energy is channelled from 
waves to protons, therefore causing an increase with time 
of the relative weight of secondary electrons and positrons 
with respect to relic electrons. 

The evolution of the spectra of particles (protons and 
electrons) and Alfven waves is obtained by solving numer- 
ically Eqs. 10, 11, and 12 with Qe{Pe,t) given by Eq. 36. 
The spectrum of the secondary electrons and positrons at 
the beginning of the acceleration period is computed from 
the Fokker-Plank equation (Eq. 10, with the source term 
given by Eq. 36) under stationary conditions and assuming 
Dpp = (e.g., Dolag & Ensslin 2000): 



Neip) 



1 



(dp) :(dp) 



Qe{p)dp. 



(55) 



Our general findings, illustrated in Fig. 2, are summa- 
rized below. 

As in the case of the reacceleration of relic primary 
electrons (Paper I), the efficiency for lepton acceleration 
decreases with increasing energy of relativistic hadrons in 
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Figure 2. Time evolution of the spectra of leptons. Tlie spectra are given at t = (tiiick solid lines), 3 X 10^^ (dashed lines), 7 X 10^^ 
(dotted lines), and 1.2 x 10^'' seconds (solid lines) after the start of the acceleration phase. Panel a): calculations are performed assuming 
rifii = 2 X 10~^cm~^ and Bp = 8% of the thermal energy density. Panel b): calculations are performed assuming n^t^ = 10~^cm~^ and 
Sp = 8% of the thermal energy density. Panel c): calculations are performed assuming nth. = 10~*cm~^ and £p = 1% of the thermal 
energy density. For all the panels the other parameters are s = 2.2, T = lO^K, B{r = 0) = 0.5/iG and d/dt{SB)'^ /Sir = 1.3 x 10~^*erg 



the ICM. As a consequence, the prominence of the bump 
of accelerated particles that appears in Fig. 2 is expected 
to decrease when the energy content in the form of hadrons 
increases. 

A pronounced feature appears in the spectrum of lep- 
ton.s. due to the reacceleration process, namely a sharp drop 
in the spectrum, followed by a flattening. The drop can be 
easily understood in terms of balance between energy losses 
of relativistic leptons and rate of ro-oncrgization due to res- 
onant interaction with waves. The flattening is simply due 
to the secondary electrons and positrons continuously gener- 
ated in the hadronic interactions of cosmic rays in the ICM. 
Once the reacceleration period becomes longer than the ac- 
celeration time-scale, the typical energy at which this fea- 
ture appears tends to decrease with time, as a consequence 
of the enhanced damping of the waves on the proton compo- 
nent. As a general comment, we point out that the presence 
of reacceleration boosts the number of leptons with 7 ~ 10* 
by 1-2 orders of magnitude with respect to the case without 
reacceleration. 

After the end of the reacceleration stage, the spectrum 
of protons remains basically unchanged. As a consequence, 
the spectrum of secondary electrons and positrons generated 
after the reacceleration is stopped is also time-independent, 
and determined only by the duration and by the efficiency 
of the reacceleration phase. 



4.5 Reacceleration of Relic and secondary 
electrons: Hybrid Models 

In Paper I we discussed at length the possibility of reacceler- 
ating relic electrons, possibly injected within the cluster vol- 
ume either at shocks, or in AGNs or ordinary galaxies. When 
these injection processes take place at redshift 2 < 0.5, elec- 
trons with Lorentz factor 7 ~ 200 — 500 may have a lifetime 
of ~ 10^ to 10^" yrs depending on whether they are injected 
in the center or in the outskirts of the cluster, respectively 
(Sarazin 1999; Brunetti 2003). It follows that, at least in the 



external parts of a cluster there may be a sufficiently high 
abundance of relic electrons to be reaccelerated by the MHD 
turbulence. 

In this Section wo describe the result of the reacceler- 
ation of both components, namely relic electrons and sec- 
ondary electrons and positrons. We name these scenarios 
Hybrid Models. 

As in the previous Section the evolution of the spectra 
of particles (protons, electrons and positrons) and Alfven 
waves is obtained by solving numerically Eqs. 10, 11, and 
12 with Qe{pe,t) given by Eq. 36. In this case, however, 
the initial spectrum of the electrons is the combination of 
the spectrum of the secondary electrons and positrons at 
the beginning of the acceleration period (Eq. 55) and of the 
spectrum of the relic-primary electrons accumulated in the 
ICM. 

In Fig. 3 we plot the spectrum of the reaccelerated elec- 
trons and positrons as obtained in the cluster center and in 
the cluster outskirts (see the caption for the numerical values 
adopted there). In Fig. 3 we assume that the strength of the 
magnetic field in the cluster volume scales according with 
flux conservation {B oc n'^l^) and that the injection power 
of Alfven waves is Pa = / Ikdk oc (see Sect. 5.1 and Pa- 
per I). The number density of the relic (primary) electrons 
and of the relativistic hadrons in the cluster volume scales 
with that of the thermal particles. Assuming that the en- 
ergy density of the relativistic protons is of the order of 10~^ 
times that of the thermal plasma, from Fig. 3 it follows that 
the reacceleration process of secondary electrons/positrons 
could be important. In particular, we find that a relevant 
contribution (from ~50 to 80%) to the spectrum of the ra- 
diating electrons and positrons with 7 ~ 10^ — 10* in the 
central regions of the cluster is provided by reaccelerated 
secondary electrons and positrons. On the other hand, the 
spectrum of electrons in the external regions is essentially 
contributed by the reaccelerated relic electrons. In the pres- 
ence of processes that make the spatial distribution of rela^ 
tivistic protons broader than the thermal gas, we can expect 
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Figure 3. Spectra of accelerated electrons and positrons in a 
cluster core (upper curves) and in the cluster outskirts (lower 
curves) at t = 7 x 10^^ s from the beginning of the accel- 
eration phase. Calculations in the cluster core are performed 
assuming n^^ = 2 X 10~^cm~^, B(r = 0) = 1.5/iG, and 
d/dt{SB f /Stt = 1.19 X 10"28gj.g g-i cm-3. In the cluster out- 
skirts we assume riff^ = 7 X 10^'''cm^'"', B(r = 0) = 0.16/iG, and 
d/dt{SB)'^/8n = 7.8 X IQ-^Ocrg s-i cm-3. = i% of the ther- 
mal energy density (with s = 2.2), T = 10*K are assumed in 
both the cluster core and outskirts. The thin solid curves axe the 
total spectra derived from the reacceleration of a population of 
relic primary electrons with fg = 10~®£th mixed with the popu- 
lation of secondary electrons and positrons. The thick solid curves 
mark the contribution of the accelerated relic electrons only. For 
comparison, the stationary spectra of the secondary electrons and 
positrons in the ICM at the beginning of the acceleration are also 
given in both cases (dotted lines). The high energy cut-off in the 
spectrum of secondary particles is due to a high energy cut-off 
fixed at 10 TeV in the spectrum of protons. 



that the contribution provided by the central denser regions 

of the cluster gets somewhat suppressed. At the same time, 
since in the cluster outskirts the Coulomb losses of the rel- 
ativistic electrons are less efficient than in the cluster cen- 
ter, the radial distribution of the number density of relic 
electrons may be broader than that assumed in Fig. 3 and 
consequently their contribution to the total spectrum of the 
reaccelerated leptons may be even larger. 



5 HYBRID MODELS: PHENOMENOLOGY OF 
NON THERMAL EMISSION FROM 
GALAXY CLUSTERS 

Here we apply the formalism described in the previous Sec- 
tions and calculate the expected non thermal omission from 
gala^xy clusters, when both relic electrons and secondary 
electrons and positrons are present during the stage of in- 
jection of turbulence and resonant reacceleration. 



5.1 Basic Assumptions 

In this Section we briefly discuss the basic assumptions 
adopted for the calculations of particle acceleration and 
non-thermal emission from galaxy clusters. The assump- 
tions are relative to the physical properties of the ICM and 
of the relativistic component, and to the modelling of the 
injection of turbulence in the cluster volume. 

We assume a Id-model (Cavaliere & Fusco-Femiano, 
1976) for the radial density profile of the thermal gas in 
the ICM, in the form 

/ r ,\ -3/3/2 
nthir) = nth{r = 0)[l + {-fj , (56) 

where rc is the core radius and /3 = 0.8. The magnetic field 
is taken in its flux conserving form: 

^W = ^(-0)(^^f^)^^^ (57) 

Here we explore the region of values B{r = 0) ~0.5-3/LtG, 
expected to reproduce the uncertainty in the value of the 
magnetic field as derived from different techniques. 

Following Fujita et al.(2003) and paper I, we assume 
that large scale fluid turbulence is injected in the ICM during 
cluster mergers and that the turbulent eddies at small scales 
radiate MHD waves due to the Lighthill mechanism (Sect. 
3.2.3). 

For simplicity we assume that the maximum injection 
scale of the turbulence, the Reynolds number and the veloc- 
ity of the turbulent eddies, which are essentially unknown 
quantities, arc independent from the location within the 
cluster volume. Under these simplified conditions, in Paper 
I we showed that the injection power in the form of Alfven 
waves scales as: 

PA{r) = [ h{r)dk = PA{r = 0)f^f^y^\ (58) 

Finally we assume that the spatial profile of the number 
density of the relic electrons and of the relativistic protons 
(at the beginning of the acceleration period) scales with that 
of the thermal matter: 

fp[e] = SthTjple], (59) 

where r)p[f,] is a parameter; a reference value, s = 2.2, is 
adopted for the injected spectrum of the cosmic ray protons. 
As a matter of fact the use of Eq.(59) limits ourselves to as- 
sume a fairly uniform distribution (non-patchy) of relativis- 
tic particles in the ICM. On the other hand, patches of rela- 
tivistic plasma separated from the thermal pool are clearly 
detected by radio and X-ray observations (e.g., Fabian et 
al. 2000; McNamaxa et al., 2000) which show the presence 
of bubbles and cavities in a few galaxy clusters. However, 
it is expected that such bubbles will expand with time and 
mix in the ICM due to the developing of instabilities in a 
time scale of the order of a few lO^yrs (Churazov et al., 
2000; Briiggen & Kaiser 2001) which is a short time scale 
with respect to the duration of the injection process of the 
bulk of cosmic rays in galaxy clusters. Thus the assumption 
of a fairly uniform mixing between thermal and relativis- 
tic plasma in galaxy clusters is justified for the aim of the 
present paper. 
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5.2 Particle acceleration in cluster cores 

In this Section we focus on the synchrotron emission ex- 
pected from the cores of galaxy clusters and derive con- 
straints on the physical properties of the ICM. We focus 
on cluster cores since in these regions the larger gas density 
makes the pp interactions more frequent and therefore the 
density of secondary particles larger. Moreover the magnetic 
field is expected to be larger in the center. These two facts 
imply a larger synchrotron emissivity in these regions. On 
the other hand, if relativistic protons have more than a few 
percent of the local thermal energy density, the damping of 
waves becomes too large and the electron acceleration gets 
suppressed, therefore reducing the synchrotron emissivity. It 
follows that the general situation may be rather complex. 

This complexity is illustrated in Fig. 4 where we plot 
the synchrotron emissivity as a function of the ratio of the 
energy densities in the relativistic protons and the thermal 
plasma at the beginning of the reacceleration stage, £p/£th, 
for different values of B and of Pa- 

As expected, for small values of £p/£th the damping 
rate is less efficient than the cascading process and the syn- 
chrotron emissivity simply scales with £p/£th (Eq. 52); at 
this stage, for the adopted rates of turbulence-injection, 
the synchrotron emission is more than one order of magni- 
tude larger than in the case without reacceleration. Increas- 
ing £p/£th, the damping of waves on protons increases and 
the acceleration efHcency decreases, so that the synchrotron 
emissivity is reduced as well. Fig. 3 also shows that when 
the ratio £p/£th is larger than ~ 10 — 20%, then the reac- 
celeration of leptons with 7 ~ 10"' is beisically stopped and 
the synchrotron emissivity approaches that expected from 
the standard secondary model. 

When this saturation effect does not occur, the syn- 
chrotron emissivity in the central region may easily exceed 
the observations. As a result, such observations can be used 
to impose contraints on the physical conditions in which the 
reacceleration of secondary particles takes place. As usual 
we refer to the case of the Coma cluster as the case in which 
the observations are richer. The conclusions that we will 
draw below should not be extended to other clusters, until 
comparable wealth of data is obtained for those clusters. 

The main observational constraints can be summarized 
as follows: 

i) Radio brightness of the core region: we use the 327 
MHz VLA profile obtained after the subtraction of point- 
like sources (Govoni et al. 2001). The allowed region of the 
brightness of the core region can be estimated by subtracting 
the contribution due to the external regions to the bright- 
ness integrated along the line of sight. 

a) Radio spectrum of the core region: the 327-1400 MHz 
spectral index map of the Coma radio halo shows a plateau 
in the central regions and a prominent radial steepening 
of the slope of the spectrum (Giovannini ct al. 1993). The 
plateau region roughly coincides with the cluster core and 
the sychrotron spectrum there is slightly flatter than a ~ 
0.9. 

Hi) Radio brightness in the frequency range 2.7-5 GHz: At 
these high frequencies the radio emission imposes strong 
constraints on all flavors of secondary models, as discussed 
by Reimer et al.(2004). 



In Fig. 4 we compare our theoretical expectations, as 
obtained for different values of the model parameters {Pa, 
B and £p/£th)-, with the constraints listed above. 

The weakest limits are clearly those obtained for low 
values of the magnetic field. For B = l/iG (left panel in 
Fig. 4), low rates of injection of waves select the region with 
values of £p/£th larger than about 10~*. On the other hand, 
as soon as the rate of injection increases above ~ 8x 10~^^erg 
cm~* (which implies a total energy budget injected per 
unit volume in the form of Alfvcn waves ^ 1.5% of the 
thermal energy density), the allowed values of £pj£th drop 
below a few 10~^. The high frequency data of the Coma 
cluster impose an upper limit £p/£th < 0.05. 

Assuming B = 3/iG (right panel in Fig. 4) it is clear that 
the observed synchrotron brightness excludes big chunks of 
the parameter space. In particular, assuming an appreciable 
injection rate of energy of Alfven waves in the cluster core 
(> 0.5% of the thermal energy) it is found that intermedi- 
ate values of the ratio £p/£th are not allowed. The region 
£p/£th > 10~^ is excluded by the high frequency points in 
the spectrum of the radio halo of the Coma cluster (arrows 
in Fig. 4). 

The most difficult observational constraint to match is 
the combination of low synchrotron brightness i) and flat 
radio spectrum ii). A relatively low value of the accelerar 
tion efficiency cannot reproduce a synchrotron spectrum as 
flat as the observed one. Therefore an efficient particle ac- 
celeration mechanism is requested to boost electrons toward 
higher energies and to flatten the emitted synchrotron spec- 
trum. In order to avoid to exceed the observed brightness, 
a relatively small injection rate of secondary electrons and 
positrons is required. More quantitatively, we find that the 
parameter space with £p/£th > is excluded for this 

strongly magnetized case. 

5.3 Integrated broad band spectrum 

In this Section we illustrate our calculations of the vol- 
ume integrated fluxes of radiation generated by reacceler- 
ated electrons and positrons through synchrotron emission 
and IC. The central gas density, nth{r = 0), the (3 parame- 
ter and the core radius r,: are chosen as the representative 
values of the Coma cluster (Bricl ct al. 1992). 

In Fig. 5 we plot our results for the synchrotron spectra 
(left panel) and the IC spectra (right panel). The data points 
refer to the radio, hard X-ray and gamma ray bands. All 
curves are obtained in the assumption that the cosmic ray 
energy density at the beginning of the reacceleration stage 
is proportional to the thermal energy density at any point. 
The values of PA{r = 0) and the ratio £p/£th are not chosen 
to obtain a best fit to the data, they are only fixed in order 
to provide a viable representation of the data. 

Some general remarks emerge from the inspection of 
Fig. 5 : 

i) The synchrotron luminosity of a Coma-like radio halo 
can be easily matched even with fp/fth ~ 10"^ - 10"^ and 
~ /LiG strengths of the central magnetic field. This energy 
requirement is more than one order of magnitude below that 
of classical secondary models (Blasi & Colafrancesco 1999; 
Dolag & EnssUn 2000). 
ii) The steepening of the integrated synchrotron spectrum 
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Figure 4. Synchrotron emissivity at 330 MHz (in arbitrary units) as a function of the ratio between cosmic rays and thermal energy 
densities at the beginning of the reacceleration stage (£pl£th) for the core region of the Coma cluster calculated after a reacceleration 
period Aracc = lO'^® s. The horizontal dotted lines give the upper and lower bounds of the allowed range of the synchrotron emissivity 

obtained from the observed brightness profile at 330 MHz (estimated from the deprojection of the observed profile given by Govoni et 
al. 2001). Arrows gives the limits on the ratio Ep/Stf^ obtained from the high frequency radio data (Reimer et al. 2004) and assuming 
that all the emission at these frequencies comes from the core region. Circles mark the points in the diagram for which the synchrotron 
spectral index is consistent with that observed (from Giovannini et al. 1993). Panel a): Calculations are carried out for B{r = 0) = 1/iG. 
The different curves refer to different assumed injection rates for Alfven waves {d/dt{6B)^ /8iv): 1.87, 1.3, 0.83, 0.64, and 0.47 x 10~^*erg 
s~^ cm~^ (from top to bottom, respectively). Panel b): Calculations are carried out for B(r = 0) = 3//G. The different curves refer 
to different assumed injection rates for Alfven waves {d/dt{5B)'^ /S-k): 1.3, 0.83, 0.47, 0.33 and 0.21 x 10" ^^erg s ^ cm ^ (from top to 
bottom, respectively). 



of the Coma radio halo can be reproduced only for cosmic 
ray energy density of the order of £p < 5% or lower. This 
is because a larger content of cosmic ray protons would de- 
crease the efficiency of the lepton acceleration and reduce 
the synchrotron bump at lower frequencies. 

iii) The observed hard X-ray spectra arc hardly achievable 
if only the effect of reaccelerated secondary particles is taken 
into account. This is mainly due to the low number of sec- 
ondaries generated if the energy density in the form of relar 
tivistic protons is the one inferred in our point i). 

If the injection of waves takes place on a spatial scale 
which is appreciably broader than that in Eq. 58, the effi- 
ciency of the reacceleration of secondary particles increases 
in the external volume and the IC omission in the outskirts 
can be enhanced, leaving the synchrotron emission almost 
unaffacted, due to the rapid decrease of the value of the mag- 
netic field with radius. However, we find that a fiux of HXR 
close to the observed one can be obtained only by assum- 
ing rather extreme conditions in the cluster outskirts (e.g., 
£t ~ £th)- Furthermore we find that in this case the strong 
back-reaction of the accelerated protons would suppress the 
acceleration of electrons within ~ lO^yrs. 

Similarly an appreciably larger IC luminosity cannot be 
produced if the spatial distribution of the cosmic ray pro- 
tons is broader than that of the thermal plasma. Indeed, 
although in this case a larger number of secondary elec- 
trons and positrons is produced in the external volume, the 



stronger back reaction of protons on the waves inhibits the 
acceleration of electrons/positrons in these regions. 

iv) Due to the relatively poor efficiency of the reaccelera- 
tion mechanism, Alfven waves in the ICM cannot accelerate 
very high energy electrons (say 7 > 10^) and thus the en- 
ergy distribution of the electrons and positrons which arc 
responsible for the emission of gamma rays though IC is not 
appreciably affected by the reacceleration scenario discussed 
here. This is the reason why in the pure reacceleration mod- 
els, with no protons and no secondary particles, we expect 
no gamma ray emission. The limit imposed on the energetic 
budget in the form of high energy protons from the EGRET 
upper limit (Reimer et el. 2004) is at the level of 20% of the 
thermal energy, obtained by assuming that the all gamma 
ray flux is generated through pion decays. For cosmic ray en- 
ergy densities below this bound (but larger than a few times 
10~^ of the thermal energy), an anti-correlation appears be- 
tween the IC HXRs and the IC gamma rays (we recall that 
the latter are only generated by the secondary electrons and 
positrons) ; a similar anti-correlation is expected between the 
IC HXRs and the gamma rays generated by the decay of 
7r°. A detailed analysis of the gamma ray emission expected 
from reacceleration models will be presented in a forthcom- 
ing paper (Brunetti et al., in prep.). 

t)j In principle the broad band non-thermal spectrum of 
the Coma radio halo can be reproduced if both relic elec- 
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Figure 5. Expected synchrotron (left panel) and IC (right panel) spectra calculated for the Coma cluster. All the calculations are 
carried out assuming that the injection rate of Alfven waves in the center of the cluster is Pa('' = 0) = 1.3 x 10~^*erg s~^ cm~^ and for 

^tacc = 10^^ s. Calculations arc carried out for both Pa oc n^j^^ (Eq. 58) (thin lines) and Pa (x ri'J^'^ (thick lines). The different models 
refer to the following cases: B{r = 0) = 2.0/iG and £p/£th = 0.12 with no relic primary electrons (dotted lines), B{r = 0) = 1.7/iG 
and £p/£tfi = 0.03 with no relic primary electrons (solid lines), and B{r = 0) = 1.1//G and Ep/Eth = 0.01 with E^/S^f^ = 5 X 10~® 
independent of the distance from the center (dashed lines). 



trons and secondary electrons/positrons are present^. If a 
few percent of the thermal energy arc stored in the relativis- 
tic protons, the synchrotron spectrum may be dominated 
by the radiation from reaccelerated secondary electrons and 
positrons, mainly in the central regions of the cluster. At the 
same time, the IC emission is dominated by reaccelerated 
relic electrons in the external regions (i.e. at 2 < r/rc < 5). 

Obviously the contribution of the reaccelerated secondary 
electrons and positrons to the integrated non-thermal spec- 
trum of galaxy clusters is expected to fall down if the energy 
of cosmic ray protons is maintained well below ~ 1% of the 
thermal energy. 

5.4 Radial profiles 

The very broad extension of the synchrotron emission from 
giant radio halos is among the properties which are diffi- 
cult to be fitted by secondary models (e.g., Brunetti 2004 
and ref. therein). Although this Section is not devoted to a 
detailed comparison between observed synchrotron profiles 
and model expectations, here we show that Alfvenic reac- 
celeration of secondary electrons and positrons in the ICM 
may generate relatively broad synchrotron profiles. 

If the reacceleration period is much longer than the 
reacceleration time-scale, the bulk of the secondary elec- 
trons and positrons injected above the momentum, p , at 
which Coulomb losses outweight the acceleration efficiency, 

t For simplicity we assume Eq. 59 to hold also in the case of relic 
electrons. 



is essentially boosted around a maximum momentum, pmax, 
at which acceleration is balanced by radiative losses. 

From Figs. 12 & 13 of Paper I one finds that under the 
assumed physical conditions the acceleration of the lower en- 
ergy electrons and positrons typically happens in the regime 
Ts » Td while the acceleration of the higher energy leptons 
happens in the opposite regime. 

Thus, assuming for simplicity a power law energy distri- 
bution of the relativistic protons Np = Kpp~" , from Eqs.(l), 
(8), (51), and (53) one finds: 



Ac nthfp 



-*) 



l/(w-Hs-2) 



AlfJ Iq 

where Ac is the constant in Eq. 1, and 



2(l-|-w-|-s) (e/c 



|3— (Lt)-t-s) 



Plow/w-pC 



s-2 (w -h s)2 - 1 
while from Eqs.(2), (8), (49), and (53) one has: 



Pmax — 



nth 



(60) 



(61) 



(62) 



where A^aA is the constant in Eq. 2, B?c-h = + and 



j4ww — _ 7r6 



1/3 



5(w - 1) 



2/3 



(63) 



The number density of the reaccelerated secondary electrons 
and positrons around pmax can be estimated by the integral 
of the number density of secondary particles injected with 
p > during the reacceleration stage. From Eqs.(38) (ne- 
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glecting for simplicity the contributions from the 7^2 and Vs 
terms) and (60), one has : 



+ ATA(g)c-Pi 
N ± ~ nth if p 



Ac nti^SpB'- 



.(64) 



On the other hand, the number density of electrons with 
p ~ Pmax in the classical secondary model can be obtained 
from Eqs.(38), (55) and (62) : 



(65) 



Thus the increase of the extension of the emitted syn- 
chrotron profile (associated to electrons with p ~ Pmaxc) 
in the secondary-reacceleration model with respect to that 
in the classical secondary model can be directly estimated 
from the ratio: 



nth 



(66) 



which, in the assumption that £p oc nth, roughly scales as 



I 1/2 3/2 

A'Ti^ /Ng± oc n^^^ - n^^^ . Eq. 66 provides a simple way 
to estimate the increase of the extension of the radial pro- 
file (at least for s ~ 2 — 2.3). However, the injection spec- 
trum of secondary electrons/positrons at lower energies (i.e., 
7 < 10'^) is not well described by a simple power law (Eq. 
38), and the spectra of protons as affected by reacceleration 
are not power laws either. Thus a more detailed calculation 
is required. In Fig. 6 we plot the ratio of the synchrotron 
emissivities as given with reaccelerated particles and in the 
context of classical secondary models (the parameters are 
as used in Fig. 5). It is clear that the models invoking reac- 
celeration of secondary particles may produce broader syn- 
chrotron emission with respect to secondary models. How- 
ever, we also notice that the presence of a cut-off in the 
synchrotron spectrum as obtained in reaccelerated models 
(which occurs at lower frequencies with increasing radius) 
produces a steepening of the synchrotron spectrum which 
balances the increase of the ratio Jgy„/ Jsyn with the dis- 
tance from the center. 

In order to reproduce the extension of the largest radio 
halos with the assumptions in Sect. 5.1, the synchrotron 
emissivity in the external regions should be amplified by a 
factor of ~ 100 with respect to that produced in classical 
secondary models (e.g., Brunetti 2004, Fig. 2); this large 
factor is not obtained in our calculations (Fig. 6). In order 
to obtain a very broad synchrotron radio emission, the radial 
distribution of the energy density of cosmic ray protons can 
be forced to be more extended than that of the thermal 
ICM. However, in this case the synchrotron profile is not 
broadened enough because at large distances the damping 
of the Alfven waves due to cosmic rays gets stronger. 



6 DISCUSSION AND CONCLUSIONS 

We presented the results of the first self-consistent calcula- 
tions of the injection and reacceleration of cosmic ray pro- 
tons and both primary and secondary leptons in the intr- 
acluster medium. This goal is reached by solving the cou- 
pled equations for the time-dependence of the protons, elec- 
trons/positrons and Alfven waves spectra. The work pre- 




Figure 6. Normalized ratio between the synchrotron emissiv- 
ity (at 330 MHz) produced by reaccelerated secondary particles, 
Jgy^, and that produced by the stationary injection of secondary 
particles, Jsym ^ ^ function of the distance from the cluster cen- 
ter (in units of the core radius) . The reacceleration models are the 
same as in Fig. 5 (with the exception of the hybrid models shown 
as dashed lines in Fig. 5) and arc shown with the same symbols. 
The synchrotron emissivity in the classical secondary model is 
computed in the assumption that B{r = 0) << 3/iG. The dot- 
dashed lines reproduce the behaviour oc n^^^^'^ and oc n^^^^ as 
expected from Eq. 66. 



sented here is the natural continuation of a previous work 
(Brunetti et al. 2004). In particular, in the present paper 
we illustrate novel results concerning the role of the injec- 
tion and further re-energization of secondary particles due 
to their resonant interaction with Alfven waves. Such sec- 
ondary particles must exist in the ICM, at least to some 
extent to be defined: the reason for such a certainty is that 
we are aware that there are cosmic ray sources in clusters 
and that the bulk of cosmic rays generated by such sources 
are diffusively trapped within the ICM (Volk et al. 1996; 
Berezinsky, Blasi & Ptuskin 1997). However, neither the 
amount of energy in the form of cosmic rays, nor their spec- 
trum or their spatial distribution within the ICM are cur- 
rently known. Their production could be related to the for- 
mation of the large scale structure of the universe or to the 
astrophysical sources within clusters (galaxies and AGNs). 

The main assumption of our calculation is the choice of 
a specific type of MHD waves, namely Alfven waves. In gen- 
eral one could think of any combination of different modes 
which behave differently in their interaction with particles 
(electrons and protons). In particular, the case of magne- 
tosonic waves, considered by Cassano & Brunetti (2005) ap- 
pears to be particularly interesting since in that case the 
problem of wave injection on small spatial scales is allevi- 
ated. In the present paper, the presence of Alfven waves 
on small scales is solved by assuming that fluid turbu- 
lence is injected on large scales and gets eventually coupled 



© 0000 RAS, MNRAS 000, 000-000 



Alfvenic reacceleration of relativistic particles in galaxy clusters in the presence of secondary electrons and positron 



with Alfven turbulence on small scales through the so-called 
Lighthill mechanism. Obviously, direct proofs that this pro- 
cess may be at work require MHD simulations with resolu- 
tion and dynamical range which are not available at present. 

The wavo-particlc coupling implies a relevant effect on 
the presence of secondary particles, mainly in two ways: first, 
the spectrum of parent cosmic rays is changed by the reac- 
celeration process (as also found in Paper I) therefore af- 
fecting the spectrum of secondary electrons and positrons. 
Second, the secondary electrons and positrons are in turn 
re-energized by the resonant interaction with Alfven waves. 
The effect of the reacceleration on the confined cosmic rays 
is shown in Fig. 1 and most notably consists of a bump at 
Lorentz factors below ~ 10^, which becomes increasingly 
more evident with time after the start of the reacceleration 
phase. At larger Lorentz factors the spectrum of cosmic rays 
remains unchanged. A similar bump shows in the spectrum 
of the total population of electrons and positrons (primaries 
plus secondaries) as plotted in Fig. 3. The pumping of energy 
into relativistic particles through the resonant interactions 
with Alfven waves at some point produces an interesting 
effect, which in Paper I we named wave-proton boiler, the 
reacceleration of electrons and positrons continues provided 
the wave damping on the proton component is not too large; 
when the energy present in the form of protons exceeds some 
limit the reacceleration of electrons gets suppressed. In this 
sense the system made of protons, electrons and waves is 
self-regulated. 

One of the most common criticisms to the so-called reac- 
celeration scenarios is that the origin of the relic electrons 
to start with is left as an open issue: we think that the 
presence of the contirmously generated secondary particles 
largely mitigates this problem, since pp interactions contin- 
uously inject these particles in the ICM. Moreover, astro- 
physical sources such as the lobes of radio galaxies and ac- 
tive galaxies are seen to pollute the ICM with a population 
of relativistic electrons, therefore it is plausible that both 
leptons of primary and secondary origin may be present in 
a cluster. The secondary particles are expected to be more 
abundant in the denser central parts of clusters, where the 
magnetic field is also larger and a correspondingly large con- 
tribution to the radio brightness is expected. In the outskirts 
of the cluster, secondary electrons/positrons are more sparse 
because of the lower gas density and a dominant primary 
component should emerge. Therefore the phenomenology of 
the non thermal activity in a cluster may reasonably be ex- 
pected to be quite complex. Contrary to the classical mod- 
els of reacceleration of relic primary electrons, the models 
discussed in this paper do predict that clusters may also 
be sources of gamma rays both due to the decay of neu- 
tral pions and to radiative processes of secondary electrons 
and positrons, whose spectrum at energies above those at 
which a reacceleration bump is generated is left basically 
unchanged. 

We applied our calculations to a Coma-like cluster in 
order to check if the main phenomenological aspects of the 
non-thermal activity may be reproduced. In doing that we 
are forced to make the simple assumption that MHD turbu- 
lence uniformly fills the cluster volume and that relativistic 
particles are efficiently mixed with the thermal pool. Clearly 
more detailed calculations are desirable, possibly making use 
of next generation numerical simulations. 



At variance with standard reacceleration models, in 
which the source of relic electrons is essentially a free pa- 
rameter, here the injection of secondary particles is self- 
regulated by the relativistic hadrons and by the reaccelera- 
tion process itself, and this increases the predictive power of 
the model. We find that the synchrotron emissivity at 330 
MHz from the core region allows us to infer useful bounds 
on the energy content in the form of relativistic protons: 
in general, as illustrated in Fig. 4, observations are hardly 
explained unless the energy density in cosmic rays at the be- 
ginning of the acceleration stage is typically less than about 
0.5% of the thermal energy density in the cluster center (at 
least when both the central brightness and the relatively flat 
spectral index, a ~ 0.7 — 0.9, of the core of Coma C are con- 
sidered). The rate of injection of energy in the form of Alfven 
waves is effectively constrained and depends on the ratio 
£p/£th- For £p/£th larger than a few times 10~* the injection 
rate is required to be of the order of ~ 10~'^*erg s^^ cm~^ 
(roughly 2% of the thermal energy during the reacceleration 
time). A larger rate implies an even lower allowed fraction of 
cosmic rays in the ICM, while an injection rate significantly 
smaller would produce a synchrotron spectrum steeper than 
that observed. Clearly the limits are more stringent when 
the magnetic field is larger. Similarly these limits are ex- 
pected to be more stringent when the real injected spectrum 
of cosmic ray protons is considerably steeper than s = 2.2 
(adopted throughout the paper) since in this case the back- 
reaction of protons on the waves is even stronger (Paper I). 
Finally, for acceleration periods longer than a few lO^yrs 
(adopted in Sect. 5) the limits are expected to be slightly 
more stringent since in this case the energy budget of pro- 
tons increases, while the limits on £p/£th would be slightly 
less stringent for shorter reacceleration stages (provided that 
they are longer than the typical acceleration time). 

We studied in detail how the volume integrated syn- 
chrotron spectrum and the radial profile of the radio emis- 
sion change when the contribution of reaccelerated sec- 
ondary electrons and positrons is taken into account. The 
dependence of the results upon the choice of the several pa- 
rameters involved is illustrated and summarized in Sections 
5.3 and 5.4. In principle given viable assumptions the radio 
emission of Coma (volume integrated spectrum and radial 
profile) could be explained in terms of reaccelerated sec- 
ondary electrons and positrons only. In particular, we have 
shown that the reacceleration of secondary leptons could 
produce synchrotron profiles which are broader than those 
expected from the standar acceleration model, although the 
broadening is still not sufficient to explain the most extended 
radio halos (at least under the assumptions in Sect. 5). On 
the other hand the observed, though controversial hard X- 
ray excess does require an additional component, that can 
plausibly be associated with primary electrons in the out- 
skirts of the cluster, where secondary particles are not abun- 
dant enough. 

The possibility to demonstrate the existence of cosmic 
ray hadrons and secondary electrons/positrons in the ICM 
is related to the detection of gamma rays due to the decay 
of the neutral pions and of circular polarization (provided 
this is not shrouded by the effect of Faraday depolarization 
in the ICM). In this respect the Hybrid Models presented 
here are qualitatively similar to the secondary models. On 
the other hand, one way of testing a scenario in which non 
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thermal radiation is generated by reaccelerated leptons of 
both primary and secondary origin is to look for a radio tail 
that should be produced through synchrotron emission of 
those secondary leptons that are not affected by resonant 
interactions with Alfvcn waves (namely with 7 > 10'"'). This 
radio flux should appear as a sort of recovery of the radio 
emissivity at high frequency, above the cutoff typical of the 
reacceleration scenarios. This flux level might be accessible 
to next generation radio telescopes, such as SKA, provided 
the radio flux at such high frequencies is not overwhelmingly 
smaller than the microwave flux at the same frequecnies. 

In this paper we focus on a specific aspect that is the 
effect of MHD turbulence on the acceleration of both pri- 
mary and secondary rclativistic particles in the ICM and 
on the related non-thermal emission from galaxy clusters. 
If future observations of radio halos will confirm the pres- 
ence of distinctive features (for example spectral steepen- 
ings) which are expected in the case of in situ acceleration 
of the emitting particles, MHD turbulence in the ICM would 
be unavoidable. In this case, detailed calculations of parti- 
cle acceleration will be important to indirectly constrain the 
presence and the properties of the MHD turbulence in the 
ICM and hopefully its connection with the magnetic field 
amplification. 

In general, turbulence could play a role in several as- 
pects of the physics of the ICM. Large-scale turbulent mo- 
tions in the ICM may provide a substantial pressure support 
to the ICM (Kulsrud et al. 1997; Rocttiger ct al. 1997; Ricker 
& Sarazin 2001), and in addition to other proposed mecha- 
nisms (e.g., Bohringer et al., 2002; Ciotti & Ostriker, 2001), 
the dissipation of turbulent energy can provide a source of 
heating to balance the cooling of cluster cores (Fujita, Mat- 
sumoto & Wada 2004) . The knowledge of the basic aspects of 
MHD turbulence in gala:xy clusters is also crucial to model 
the transport of heat and metals in the ICM (Cho et al., 
2003; Voigt & Fabian 2004). 

Future experiments, such as ASTR0-E2 (and NEXT), 
would hopefully constrain the energy budget associated to 
the turbulent eddies in the ICM by looking at the profile of 
the FeK-lines (and other) in the X-ray spectrum of galaxy 
clusters (Sunyaev et al., 2003). 
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